Qubit state detection using the quantum Duffing oscillator 
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We introduce a detection scheme for the state of a qubit, which is based on resonant few-photon 
transitions in a driven nonlinear resonator. The latter is parametrically coupled to the qubit and 
is used as its detector. Close to the fundamental resonator frequency, the nonlinear resonator 
shows sharp resonant few-photon transitions. Depending on the qubit state, these few-photon 
resonances are shifted to different driving frequencies. We show that this detection scheme offers 
the advantage of small back action, a large discrimination power with an enhanced read-out fidelity, 
and a sufficiently large measurement efficiency. A realization of this scheme in the form of a persistent 
current qubit inductively coupled to a driven SQUID detector in its nonlinear regime is discussed. 

PACS numbers: 03.65.Yz, 42.50.Dv, 85.25.Cp, 42.50.Pq 



I. INTRODUCTION 

The efficient and reliable detection of the quantum me- 
chanical state of a nanoscale system is a key component 
of all present designs of quantum circuitsii One nonde- 
structive readout scheme currently in use for the impor- 
tant class of superconducting flux qubits is based on a 
heterodyne detection of the dynamic response of a dc su- 
perconducting quantum interference device (dc-SQUID) 
detector which is inductively coupled to the qubitiS~ 
Thereby, the dc-SQUID is operated in its linear regime as 
a shunted variable inductor in a resonant circuit. In this 
set-up, its resonance frequency depends on the magnetic 
flux generated by the qubit being in the ground or excited 
state. Hence, measuring the impedance of the resonant 
circuit as a function of an externally applied bias cur- 
rent yields two characteristic Lorentzian resonances at 
two different resonance frequencies, which depend on the 
two qubit states. This detection scheme, hence, allows us 
to infer the state of the qubit from the resonant response 
of the detector in the nanocircuit. In order that a reliable 
discrimination of the two qubit states becomes possible 
in this continuous type of readout design, the probabil- 
ity distributions for the readout values have to be only 
weakly overlapping. Due to thermal and quantum fluc- 
tuations, the readout naturally is a random process;^ and 
the noise properties of the nanocircuit around the detec- 
tor resonances determine the discrimination power of the 
set-up. 

An alternative readout scheme is the Josephson bi- 
furcation amplifierj^i^ It is based on a classical driven 
nonlinear resonator and exploits the classical bifurcation 
point of the dynamically induced bistability with a small- 
and a large-oscillation statei^ The response (or output) 
of the nonlinear resonator around the bifurcation point 
is very sensitive to small changes in the circuit parame- 
ters. This is an ideal prerequisite for a sensitive detec- 
tor. Depending on the state of the qubit to be sensed, 
the resonator bifurcation point is shifted to a different 
frequency, allowing for large discrimination powers be- 
tween the large- and small-oscillation detector state of 



up to 98%fi Nevertheless, since the detector is a classical 
macroscopic device, it introduces considerable dephasing 
and relaxation to the qubit state, yielding a reduced con- 
trast of the qubit Rabi oscillations of less than 90% This 
implies that the thermal noise properties of the nonlin- 
ear detector (together with semiclassical corrections due 
to quantum fluctuations) around the classical bifurca- 
tion point determine the discrimination power between 
the two states close to the classical bifurcation point4"— 
Hence, it would be desirable to combine the advantage of 
a large discrimination power of a nonlinear detector with 
the reduced noise sensitivity of a nanocircuit operated 
close to the quantum regime. 

In this paper, we introduce a combination of both 
strategies and propose a nonlinear detector scheme in the 
form of a nonlinear resonator with an amplitude modu- 
lated drive in its few-photon deep quantum regime. In 
particular, in this regime, we shall exploit sharp multi- 
photon resonances in the nonlinear resonator^"— which 
are induced by the external driving field close to the fun- 
damental resonator frequency. They can be used for the 
detection of the states of the qubit and offer the advan- 
tage of being rather sharp and externally tunable by vary- 
ing the parameters of the external drive. The concept is 
an extension of the case of a linear resonator, where the 
fundamental resonance frequency is shifted depending on 
the qubit state. However, the multiphoton resonances in 
the nonlinear detector close to the detector's fundamental 
frequency show very small line widths. The width of the 
A'^-photon resonance is determined by the corresponding 
A'^-photon Rabi frequency, which decreases with increas- 
ing photon number. The sharp resonance lines, in turn, 
offer the advantage that only a few measurement cycles 
are necessary to ensure a large discrimination power. To 
understand the back action of the nonlinear multiphoton 
detector on the qubit state, we determine the relaxation 
rate of the qubit due to the coupling to the driven dissipa- 
tive nonlinear oscillator around a multiphoton resonance. 
Notably, the back action of the resonator on the qubit is 
sufficiently weak, yielding to a good qubit-state measure- 
ment fidelity. Furthermore, we show that the discrimina- 
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tion power of the set-up is rather large and beyond 98% 
for our choice of realistic parameters of a flux qubit cir- 
cuit. In fact, it gives rise to an enhanced measurement 
fidelity as compared to the linear parametric oscillator. 
Furthermore, we show that the nonlinear multiphoton 
detector does not have a worse measurement efficiency 
as compared to the linear detector scheme. We deter- 
mine the measurement efficiency of the set-up via the 
ratio of the time it takes to collect enough information 
on the qubit state (measurement time) and the relax- 
ation time. It turns out that the measurement efficiency 
does not considerably decrease as compared to the linear 
case. Hence, the detection scheme indeed has the ad- 
vantage of an overall reduced back action in combination 
with an enhanced discrimination power, together with a 
sufficiently large measurement efficiency. 

An experimental realization of a driven nonlinear res- 
onator in its few-photon quantum regime is in principle 
possible with present set-ups and technology. In a recent 
experiment?^ a nanoscale superconducting microwave 
resonator has been driven to its nonlinear regime by fast 
frequency-chirped voltage pulses. At low enough tem- 
perature, the regime of quantum noise has been reached. 
In this experiment, the applied driving strength has been 
rather large, which corresponds to a large photon number 
transferred to the resonator. No particular few-photon 
resonances have been revealed and the nonlinear response 
is similar to previous schemes on classical bifurcation 
detectors using a time-dependent driving frequencyJ^. 
However, the route to the few-photon regime seems to 
be clear. 

The paper is organized as follows. In Sec. [TTl we 
start from a typical experimental setup for the flux qubit 
and its SQUID detector and we derive the Hamiltonian 
model. This serves to motivate an experimental real- 
ization of our proposed detection scheme. Moreover, we 
discuss the regime of validity of the model. In accordance 
with the approximation made in Sec. |TT1 we continue the 
study of the coherent dynamics in Sec. [Hi] in the rotating 
wave approximation. Dissipative coupling to the environ- 
ment is included on the level of a Born-Markov master 
equation in the rotating frame in Sec. IIVI In Sec. |Vl 
we analyze the multiphoton transitions in the nonlinear 
response of the Duffing oscillator and show that their 
resonance frequency depends on the qubit state. Then, 
in Sec. I VII we determine the back action of the driven 
dissipative detector on the qubit dynamics by analyzing 
the population difference of the qubit states at the mul- 
tiphoton transitions in the detector. Furthermore, we 
determine the measurement efficiency in Sec. I VIII 



II. MODEL 

In order to relate the theoretical approach in the fol- 
lowing to realistic devices, we start by deriving the model 
from standard set-ups already realized in experiments. 
For this, we use a typical architecture of a persistent 



current qubit which is inductively coupled to a driven 
SQUID. 



A. Persistent current qubit 

We consider the experimental set-up used in Ref. [l3 
for the qubit, consisting of a superconducting loop inter- 
rupted by three Josephson junctions, two of which have 
equal Josephson energies, while the coupling energy of 
the third is smaller, in order to yield a double-well po- 
tential configuration. In this low-inductance circuit, the 
flux through the loop remains close to the externally ap- 
plied value $qb. When the latter is close to (n -I- l/2)$0; 
where n is an integer and $o is the flux quantum, the 
device is described by the Hamiltonian in terms of the 
Pauli matrices ax,z as (ft = 1) 

i^qb = I cr^ - yo-^, (1) 

with the two eigenstates jt) and jj,) of cr^ correspond- 
ing to the two persistent current states ±/p. The min- 
imal energy level splitting A and the current Ip are 
determined by the charging and Josephson energies of 
the Josephson junctions. The asymmetry is given by 
e = 2/p(<I'qb — $o/2). In the energy eigenbasis, the Hamil- 
tonian follows as i/qb = UJqhTz/2, with Wqb = 

and Tr = CTr cos 6 — ax sin 9 is the corresponding Pauli 
matrix with tan^ = A/e. The detection of the qubit 
state essentially involves the measurement of the mag- 
netic flux produced by the persistent current states. To 
this end, one can use the driven SQUID as a sensitive 
magnetometer^ operating in its nonlinear region. Below, 
we will restrict to the few-photon deep quantum regime. 



B. Driven SQUID as a nonlinear quantum detector 

We consider the standard setup of a dc-SQUID formed 
by two Josephson junctions in a superconducting loop, 
but subject to a time-dependent external bias currenti^ 
Moreover, we assume a negligible ring inductance Ln 
of the SQUID (low- inductance approximation) In this 
conflguration, the superconducting phase differences at 
each junction, xi and X2, play the role of dynamical vari- 
ables with a constraint given by the flux quantization, 
i.e., xi - X2 = -^sq/tpo = -27r(y9cx, where <I>sq is the ex- 
ternal magnetic flux piercing the superconducting loop 
and ipo — ^o/Stt. Note that within the low-inductance 
approximation, Li^/oc with the critical current /qc 

of the SQUID. Thus, the system is described by the gen- 
eralized coordinate x+ = (Xi + X2)/2, with the effective 
LagrangianiS 

Lsq{x+,X+,t) = Cq x+ + Ej cos {iripcx) cos(x+) 

-(^o/b(i)x+, (2) 
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where we have assumed a symmetric loop, with Ej = 
ipoloc as the Josephson energy, and Cq as the capaci- 
tance of each junction. Moreover, we include a time- 
periodic ac current Ih{t) = locos^ujc^t) with frequency 
and amplitude Iq injected "into" the loop. The 
above Lagrangian describes an effective superconduct- 
ing loop (with a negligible ring inductance) with a sin- 
gle Josephson junctioi*^ with a tunable Josephson energy 
-Ej cos(7r(/3cx); critical current Ic = 2/oc | cos tti/Sox | , cross- 
junction phase difference x+, and capacitance C = 2Co. 
In order to tune the resonance frequency, the SQUID is 
shunted^ with a capacitance Cs ^ C. Next, we shall 
establish the optimal working point of the qubit-detector 
system, where the dissipative influence entering via the 
detector is minimal. 
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FIG. 1. (Color online) Dependence of the dimensionless ratios 
a/fl and g/Q on the external flux <^cx in the SQUID. The 
parameters of the SQUID are chosen as Cs — 7.65 pF, Ico ~ 
200 nA, Ip = 300 nA, and M = 40 pH.^^ 



1. Qubit-detector interaction 

The qubit and the SQUID arc coupled by means of 
their mutual inductance M 4— Thereby, the SQUID in- 
duces the flux MIq in the qubit loop, where Iq is the 
circulating current in the SQUID. The latter can be 
determined by using current conservation in the loop 
and the Josephson relations for the two junctions in 
the SQUID. For the symmetric SQUID^ it follows that 
Ic){t) = /co sin(7r(/3cx) cos(x+(i)). Thus, the total mag- 
netic flux in the qubit is affected by its coupling with 
the SQUID, and it is composed of the external flux and 
the SQUID-generated contribution, i.e., 'I'qb — > $qb + 
MJo(t). This implies that the energy bias of the qubit 
acquires a contribution that depends on the circulating 
current in the SQUID, leading to the effective asymmetry 
Cf = e($qb) + Piloit)), where /3(/o(0) = 2MJj,Jo(t). 

Therefore, two sources of noise can affect the qubit 
dynamics, i.e., the fluctuations from the external flux 
<I>qb and from the bias current Ih{t) in the SQUID 
which is related to x+ by the Josephson equation /b(t) = 
/co sin(x+(i)). By tuning the bias current to the crit- 
ical value characterized by {d/3/dI\^)j^^=j^ = 0, the 
influence from current fluctuations in the SQUID can be 
minimized^i and the optimal working point is reached. 
For a nonsymmetric SQUID, the lowest-order contribu- 
tion is linear in /h,— i^i while in the symmetric case this 
lowest-order contribution vanishes, which implies that 
around the optimal working point the phase x+ is very 
small, x+ ~ 0. In the following, we consider a setup close 
to the optimal point, where we can expand the expression 
for Iq up to second order in x+, yielding the interaction 
term 

-ffqb-sq = .gx+CTz, (3) 
with the coupling constant g — 2/p/ocM sin(7riy9cx)- 



2. SQUID modelled as a Duffing oscillator 

As we operate the detector in its nonlinear 
regime, we expand the potential term V{x+) = 
-Ej cos (7r(^ox) cos(x+) ^ Vq + mVi^x\/'^ - OiX% in Eq. 
(PI around the optimal point up to fourth order in x+^ 

where m ~ V'qC's is the effective mass. Q = (Ic/ipoCs)^^'^ 
the corresponding frequency, and a = the 
strength of the nonlinearity. We switch to a description 
in terms of the creation and annihilation operators a and 
, defined by x+ = Xo(a + a^) with xo = \/l/{2mn) 
the zero-point energy of the phase x+ ■ Adding the time- 
dependent driving term yields us to the Hamiltonian of 
the driven SQUID described by the quantum Duffing os- 
cillator model 

i7,q ^naU- ^{a + at)4 + f{a + at) cos {uj^J) , (4) 

with nonlinearity and driving strength given by a = 
SlcifioXoy ^nd / = IqipqXo, respectively. Similarly, the in- 
teraction Hamiltonian in terms of ladder operators reads 
as 

ffqb-sq = |(a -f a'l')^(T2, (5) 

with g = 2gxl- 

Notice that g and a depend on the external flux tpcx, 
i.e., they are tunable in a limited regime with respect 
to the desired oscillator frequency il, where the cou- 
pling term is considered as a perturbation to the SQUID 
{g < a), in order to keep the dynamics of the oscillator 
to dominate. The dependence of the dimensionless ratios 
a/fl and g/fi is shown in Fig. [1] We restrict to param- 
eters of the external magnetic flux in the SQUID loop, 
which generate a weak nonlinearity and a weak qubit- 
detector coupling strength, {q;,;;} ^ fi, i.e., for ipex ~ 0. 
A typical dependence of both parameters for typical ex- 
perimental parameters is shown in Fig. [1] Both cases of 
g/fl > a/il and g/fl < a/fl can be achieved. For our 
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purpose of a qubit-detector setup, the qubit-resonator 
couphng typically will be required as small enough in 
order to ensure a minimal back action. On the other 
hand, the qubit-detector coupling should be large enough 
so that an efficient detection of the qubit state becomes 
possible. As is shown in Fig. [T] and will be quantita- 
tively discussed in the sequel of this paper, this can in- 
deed be achieved for realistic parameters. Moreover, the 
choice of the parameter regime also justifies us to re- 
strict the influence of the resonator coupling on the ef- 
fective qubit bias to lowest order in x+ only. Eventu- 
ally, the total system is described by the Hamiltonian 

H{t) = TJqb -I- i?qb-sq + -f^sq(i)- 



III. COHERENT DYNAMICS AND 
ROTATING- WAVE APPROXIMATION 

Before we address the dynamics of the detection 
scheme based on the nonlinear response of the Duffing 
oscillator to the applied periodic driving in the station- 
ary regime, we discuss the coherent dynamics generated 
by H(t), which is periodic in time. 

Here, we are interested in exploiting few-photon tran- 
sitions in the detector around the fundamental detector 
frequency fl. Hence, higher harmonics have a small am- 
plitude and can effectively be neglected. Furthermore, 
we focus on the regime of weak nonlinearity, weak driv- 
ing, and weak qubit-detector coupling as characterized 
by {a, f,g} ^ fi. The proposed mechanism of detection 
is most conveniently discussed in the simplest case, when 
the dynamics occurs close to the fundamental oscillator 
resonance Wox ^ fl ^ a;qb/2. Then, the rotating- wave ap- 
proximation (RWA) can be invoked in order to obtain a 
simple interpretation in terms of few-photon transitions. 
In passing, we note that we have also performed a com- 
plete analysis in terms of full Floquct theory, thereby 
avoiding the RWA. For all cases shown below, both ap- 
proaches yield coinciding results. 

We switch to the rotating reference frame by the trans- 
formation R{t) = exp{i{a'^a + Tz)woxO- Then, the RWA 
eliminates the fast oscillating terms from the transformed 
Hamiltonian V. = R{t)H{t)R^t) - iR{t)R''{t) and the 
time-independent Schrodingcr equation in the rotating 
frame Hlipa) = ealfa) follows, with the RWA Hamilto- 
nian given by 



'R — Rqh + ^qb-sq + ^sq , (6) 

with 

^qb = ^^l^qbTz, 

'Hqb-sq = g cosO ar, + ^ smO (a^^T~ + o^t'^) , 
Hsq = 5VLa}a — —a)aaa^ + ■^^('^ + '^'^) ■ 



The detuning frequencies follow as Jfi = — Wox and 
(Jwqb = Wqb — Zwex, and = [Tx^iTy) /2. The quasiencr- 
gies Ea and the RWA eigenstates \tfa) result from a 
straightforward numerical diagonalization of H. In the 
static frame, an orthogonal (at equal times) set {|<^a(i))} 
of approximated solution of the Schrodinger equation fol- 
lows as 

\^^{t)) ~ e— °*|</.„(t)) = e— °*e-'("'"+^^)'^"-*|^„). (7) 

Here, the quasienergy states \(j)a{t)) = 
g-t(ata+T-)woxt|(^^^ are time periodic with period 
and form a complete basis that will be used 
below for the description of the dissipative dynamics. 
We note that an analytic expression for the multi-photon 
resonances would follow from a Van-Vleck pcrturbative 
approach in a similar manner as for the pure quantum 
Duffing oscillator .ii^ii^ However, the resulting expression 
will be cumbersome and not further illuminating for 
the present purpose. We note, furthermore, that the 
qubit-detector interaction occurs via a parametric 
coupling g cos 9 aTz, and via a two-photon coupling 
5sin6'(at2r- +aV+)/2. 



IV. DISSIPATIVE DYNAMICS 

The electronic nanocircuit is embedded in a dissipa- 
tive environment. In particular, the SQUID is shunted 
with an Ohmic resistor, which yields dissipative fluctu- 
ations ^(i).— We focus to the case of an undcrdamped 
SQUID, where the shunt resistance is large,^^'^ and use 
the standard harmonic bath in order to model the fluctu- 
ations, which are rooted in current fluctuations and can 
be encoded in the Ohmic spectral density J[uj) ~ ^lo^ 
They couple to the resonator's dipole operator, i.e., 
= X+^{t). We note that, in the same way, the direct 
coupling of the qubit to the electromagnetic fluctuations 
could be included. However, we have checked^^ that for 
a related set-up of a flux qubit coupled to a harmonic os- 
cillator, such a direct dissipation of the qubit yields only 
minor quantitative corrections, which should be included 
in a quantitative description of an experiment,— but do 
not add qualitatively new physics. 

The time evolution of the reduced density operator 
is described in terms of a standard Markovian mas- 
ter equation projected onto the basis of the quasienergy 
states 

Qati{^) = -«(£a - Sfi)Qaii ^ ^ Ca/S ,a' p' Qa' p' {t) , (8) 

a'l3' 

where gap{t) = {(j)a{t)\p{t)\(j)p (t)) . The dissipative tran- 



sition rates are given hy^^^ 

n 

~Saa' ^ Na"l3'^-nXp'a",-nXa"l3,n 
a" ,n 

~^l3fi' ^ Npiia' ,-nXal3" ,-nXP"a' ,n, 

(9) 

with n ^ li and Xa(3,n being the Fourier components ac- 
cording to (0a(<)IX+l</'/3(O) = I]n«^Xp{-iWox?^i}Xa/3,n- 

Furthermore, we have used the Planck numbers 
Nai3,n = N{ea ~ s fs + nuj^^), where N^e) = 
7e [coth(e/2r) - 1 + e(-e)] with ks = 1, tempera- 
ture T and 0(x) being the Heaviside function. Since 
within the rotating wave approximation, \(l)a{t)) w 
g-i(a a-i-T^)wext|(^^^^ |;j^g only non-zero Fourier compo- 
nents are Xap.i = Xo(</5a|a|</5/})/V2 , and Xap~i = 
Xo{'^a\0''^Wp)/V^ and the master equation ([8|) consid- 
erably simplifies as it involves only single step transi- 
tions, i.e., one-photon emission (for ti = — 1) into and 
absorption (for n ~ +1) processes from the bath. We 
note that neglecting also the quasienergy dependence of 
the Planck numbers would yield the well-known Lindblad 
master equation. 

In order to measure the dynamic response of the res- 
onator to the external drive at asymptotically lon g tim es, 
a heterodyne detection scheme such as in Ref. [2J can 
be used)^ where the coupled qubit-oscillator system ap- 
proaches the steady state = Q{t — > oo). We determine 
the stationary solution characterized by q{oo) = nu- 
merically. For this, we solve the corresponding eigenvalue 
problem and follows as eigenvector to the eigenvalue 
zero. With this, we compute the nonlinear response of 
the detector, characterized by the mean value {x+)oo{t) 
at asymptotic times. As we restrict the discussion to 
the regime close to the first harmonic (small detuning), 
higher harmonics can be neglected and we immediately 
obtain 

(x+>oo(t) = tr(e-x+) 

= Y.p'S0mt)\x+\Mt)) 

- 2^eafi\Xl3a, + ie +Xl3a,-ie )■ 

a/3 

(10) 

As the system is driven with frequency Wex, {x+)oa{i) 
also oscillates with time. Its amplitude is given by 

A^^Qal3{Xl}a, + l+Xlia-l)- (H) 

Correspondingly, we evaluate the population difference 
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FIG. 2. (Color online) (a) Amplitude A of the nonlinear re- 
sponse of the decoupled quantum Duffing detector (g = 0) as a 
function of the external driving frequency Wex. (b) The corre- 
sponding quasienergy spectrum Sa- The labels A'^ denote the 
corresponding A''-photon (anti-)resonance. The parameters 
area = O.Olfi, / = 0.0060, T = 0.0060, and 7 = 1.6xl0"''O. 

{<^z)oo(t) of the qubit states and obtain 
(ff.)oo(t) =tr(^,°°a,) 

= sin0^ ^,-^(r+„e2--* + V^e^^^-*) 
+ coseY,Q^f,T^p^, (12) 

all 

where r^^ = {(pa\Tz\ip0) and r^^ = {(pa\T^\(pi3). The 
population difference oscillates, with a maximal value 
given by 

Poo = cos 6* ^ Qa/3 + siu ^ QapT^a ■ (13) 

V. DETECTOR'S DYNAMICS 

A. No coupling between detector and qubit 

Before turning to the quantum detection scheme, we 
discuss the dynamical properties of the isolated detector, 
which is the quantum Duffing oscillator. A key property 
is its nonlinearity which generates multiphoton transi- 
tions at frequencies Wcx close to the fundamental fre- 
quency CI. In order to see this, one can consider first the 
undriven nonlinear oscillator with / = and identify de- 
generate states, such as \n) and \N—n) (for N > n), when 
Sn = a{N + l)/2M^^^ For finite driving / > 0, the de- 
generacy is lifted and avoided quasienergy level crossings 
form, which is a signature of discrete multiphoton transi- 
tions in the detector. As a consequence, the amplitude A 
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of the nonlinear response signal exhibits peaks and dips, 
which depend on whether a large or a small oscillation 
state is predominantly populated.—*^ The formation of 
peaks and dips goes along with jumps in the phase of the 
oscillation, leading to oscillations in or out of phase with 
the driving. A typical example of the nonlinear response 
of the quantum Duffing oscillator in the deep quantum 
regime containing few-photon (anti-)rcsonances is shown 
in Fig. HJa) (decoupled from the qubit), together with 
the corresponding quasienergy spectrum [Fig. [2]^b)]. We 
show the multiphoton resonances up to a photon num- 
ber = 5. The resonances get sharper for increasing 
photon number, since their widths are determined by 
the Rabi frequency, which is given by the minimal split- 
ting at the corresponding avoided quasienergy level cross- 
ing. Performing a perturbative treatment with respect to 
the driving strength /, one can get the minimal energy 
splitting at the avoided quasienergy level crossing (0, N) 
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Because the nonlinearity a is typically fixed by the design 
of the SQUID, the Rabi frequency can be easily tuned by 
tuning the driving strength /. 



B. Detector response for weak coupling to the 
qubit 

Next, we consider a finite coupling of the detector to 
the qubit whose state is to be sensed, i.e., g ^ 0- The 
coupling inevitably induces relaxation and decoherence 
in the qubit, characterized by the relaxation and dephas- 
ing rate, F and F^, respectively. Typically, the detec- 
tor couples only weakly to the system, i.e., g ^ Wqb. 
Then, the associated relaxation and dephasing times (Ti 
and T2, respectively) are still much larger than the corre- 
sponding relaxation time scale for the detector given by 
1/7. In passing, we note that the corresponding relax- 
ation time around a resonant multiphoton transition (in 
the underdamped case) has been shown in Refs. [l^ and 
H to be comparable to 7. Moreover, we bias the qubit 
with a large asymmetry, e ^ A in order to "gauge" the 
detector response. 

For a rough evaluation of the order of magnitude of 
the involved time scales, we may neglect the nonlinear- 
ity of the detector (a ~ 0) for the moment and es- 
timate the effective relaxation rate for the qubit cou- 
pled to an Ohmically damped harmonic oscillator)^ This 
model can be mapped to a qubit coupled to a struc- 
tured harmonic environment with an effective (dimen- 
sionless) coupling constant KcS = 87.9^ /Jl^. For the re- 
alistic parameters used in Fig. [5] and g = 0.0012f2, we 
find that Kqs ~ 10"^", giving rise to an estimated re- 
laxation rate^^i^ Fharm - (7r/2)sin^(6')Kofre ~ IQ'^^n 
(evaluated at low temperature). Hence, this illustrates 
that we can easily achieve the situation where Fharm *C 7 
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FIG. 3. (Color online) Nonlinear response A of the detector 
as a function of the external driving frequency ojcx in the 
presence of a finite coupling g = 0.0012 Q to the qubit (black 
solid line). The blue dashed line indicates the response of the 
isolated detector. The parameters are the same as in Fig. [2] 
and e = 2.2il and A = O.OSSl, in correspondence to realistic 
experimental parameters^. 



required for this detection scheme. Then, for a waiting 
time (after which we start the measurement) much longer 
than the relaxation time 7"^ of the nonlinear oscillator, 
but still smaller than F~^, the oscillator is able to reli- 
ably detect the qubit state. In fact, under these condi- 
tions, the state of the qubit, apart from the inevitable 
dephasing, remains unaffected in a time window before 
it reaches its global stationary state and an effective shift 
of the oscillator's eigenfrequency arises due to the para- 
metric coupling term ^ g cos a^a in Eq. Treating 
the qubit-detector interaction term in Eq. ^ perturba- 
tively to lowest order in g, the eigenfrequency shift follows 
straightforwardly as 

+ g (cr^) . 

Thus, the nonlinear response is shifted by —g (+.g) if the 
qubit is prepared in the state cr^ — ~1 (ctz = 1). This 
is illustrated in Fig. [3l in which we show the nonlinear 
response of the resonator for the uncoupled (blue dashed 
line) and the coupled (black solid line) case. For a fixed 
value of g, the shift between the two cases of the opposite 
qubit states is given by the frequency gap iJwox — 2 g. 
Figure 131(a) shows the nonlinear response of the detector 
for the two cases when the qubit is prepared in one of its 
eigenstates: \f) (orange solid line) and \l) (black dashed 
line). 

An important feature of a detection scheme is that it 
is efficient in discriminating the states to be detected. 
This can be quantified by the discrimination power of 
the detector, which can be defined for our case as 



A 



It) 



(15) 



The result for D{ujcx) is shown in Fig. 2] (b). The dis- 
crimination power shows a rich structure of local max- 
ima and minima, which indicates that it can be tuned 
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FIG. 4. (Color online) (a) Nonlinear response A of the de- 
tector coupled to the qubit prepared in its ground state |^) 
(orange solid line) and in its excited state |t) (black dashed 
line) for the same parameters as in Fig. |3l The quadratic 
qubit-detector coupling induces a global frequency shift of 
the response by Sujcx = 2g. (b) Discrimination power D{iUcx) 
of the detector coupled to the qubit for the same parameters 
as in a). 



directly by tuning the driving frequency. It is moreover 
important to realize that the discrimination power can 
be optimized by tuning g. In the optimized case, a lo- 
cal maximum of the multiphoton resonance for one qubit 
state can be made to coincide with a local minimum of 
the response for the opposite qubit state yielding to a 
maximal discrimination power. An example where the 
discrimination power has been optimized with respect to 
the three-photon resonance is shown in Fig.|3](b). 



VI. BACK ACTION IN THE QUBIT 

Another important prerequisite for a useful detection 
scheme is that the coupling of the qubit to the detec- 
tor around a multiphoton resonance does not generate a 
destructive back action on the qubit dynamics. In this 
section, we show that the back action in this design is 
surprisingly small for a realistic choice of parameters. 

The back action of the detector on the qubit arises 
in the form of two contributions from the coupling. 
First, this coupling has a parametric component "Hi = 
g cos OriTz, which commutes with the Hamiltonian. 
Thus, in the presence of a coupling of the oscillator to 
the bath, this term only produces dephasing and no re- 
laxation, as it is, for instance, required for a quantum 
non-demolition measurement. This part guarantees an 
efficient detection of the qubit state. The second compo- 
nent T-L2 — gsn\6{a^^T~ + a^r+)/2 in the coupling term 
yields transitions in the qubit when two-photon processes 
are induced in the detector by the external driving and / or 
by dissipative transition. Since, at low temperature, dis- 
sipation is dominated by photon leaking and the driving 
is very weak, the decay rate of the qubit from its excited 



state to its ground state accompanied by the emission 
of two oscillators photons, largely exceeds the excitation 
rate from the ground state to the excited state accompa- 
nied by the absorption of two photons originally coming 
from the bath or the driving. On the other hand, when 
the effective oscillator frequency is close to a multipho- 
ton resonance, photon absorption in the coupled system 
is enhanced and thus the asymptotic qubit population 
might be reduced. 

Thus, for a large asymmetry |e| ^ A, peaks and dips 
in the qubit population difference Poo are expected when 
multiphoton transitions in the detector are induced. This 
is what is shown in Fig. [SJa), where Poo is shown for sev- 
eral values of /. For an easier orientation, we show in ad- 
dition the corresponding stationary nonlinear response of 
the detector in Fig. [5l[b). For increasing driving, the de- 
viation from the expected value Poo = — 1 becomes more 
pronounced for larger photon numbers N and larger driv- 
ing /. The reason is that, for increasing driving, a larger 
Rabi frequency for the corresponding transition results 
[see Eq. P^ ]. From Fig. [Sj it follows that when the 
qubit is prepared in its ground state | (we consider 
e S> A) the back action is very small. The impact is less 
than 2% for the considered realistic parameters, yielding 
to a readout contrast of more than 98%. This has to be 
compared with presently achievable readout contrasts of 
less than 90% which results from an architecture with 
a classical Josephson bifurcation amplifier. In passing, 
we note that the detector response can also be calcu- 
lated from the stationary solution of the master equation 
(|5]), but for the parameters considered here (in particular 
because of the large qubit bias), this coincides with the 
shifted one. 

Moreover, we note that the components T-Li and 'H2 
can be tuned by e and A. Therefore, gsmO can in prin- 
ciple be eliminated by setting A = 0, which would im- 
ply that the measurement scheme keeps the state of the 
qubit without any relaxation but only pure dephasing 
(ideal quantum nondemolition measurement). However, 
turning off the splitting implies a major change in the 
experimental design of the sample, since this parameter 
is determined by the Josephson energy in the junctions 
of the superconducting fiux qubit and, thus, may not be 
easy to be realized. 

The back action of the detector on the qubit should be 
small not only when the qubit is in its ground state but 
also when it is in its excited state. We therefore address 
next the relaxation rate of the qubit. Energy relaxation 
in the qubit induced by the measurement process will be 
proportional to the fluctuations of the square of the phase 
operator x-i- induced by the detector's environment 
This relaxation process is characterized by the transition 
rateii^ 

Tc^f sin^eS^2^[~u^,,l (16) 
which has been computed perturbatively to lowest order 
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FIG. 5. (Color online) (a) Asymptotic population diflerence 
Poo of the qubit states, and (b) the corresponding detector 
response yl as a function of the external frequency cjox for the 
same parameters as in Fig. [3] 



in g. Here, 



5^2 [uj] 



47r 



dt 



+ 00 

dre-({x^(r + t),xi(t)}) 



(17) 

is the symmetrized power spectrum of x+ averaged over 
the period of the external driving (see Appendix for de- 
tails), with {, } indicating the anticommutator. The fact 
that information on the qubit state is acquired in the 
detector via the same channel by which dissipation is 
introduced is nicely reflected in the expression of the re- 
laxation rate in Eq. ([T5)) . In Fig. Hl^a), the relaxation 
rate F is shown for a large negative asymmetry in the 
qubit. The relaxation rate is strongly peaked around the 
multiplioton transitions. There, the noise from the de- 
tector absorbs more energy from the qubit around the 
multiphoton transition (0, N) since the parametric com- 
ponent "Hi of the coupling becomes negligible, leading to 
a dominant relaxation process induced by 7^2- 

We emphasize that although the relaxation is maxi- 
mally enhanced at a multiphoton resonance, the absolute 
value of F is still very small in comparison to the damp- 
ing constant, e.g., F/7 ^ 10~^. Thus, we can infer the 
qubit state with sufficient precision by operating the de- 
tector in its steady state regime as it has been assumed 
in Section IVBl 



time of the detector but also the time it takes to acquire 
sufficient information to infer the qubit state (the mea- 
surement time Tineas)- Hence, for a good measurement 
fidelity, Tmoas should be smaller than the characteristic 
time F-i given by Eq. (HH), or, F^eas/r > 1. 

The measurement time can be formalizedii>^i^ as the 
ratio of the symmetrized power spectrum S^^ of the 
phase operator x+ (evaluated at zero frequency) and the 
square of the difference between the two expectation val- 
ues of x+ when the qubit is in the two opposite states, 
i.e., with Eq. p3)) . 



x+ 



[Diu 



(18) 



The result for Tn 
Fig.[S]b) for the parameter set used above, for which the 
discrimination power D(u!cx) around the 3-photon reso- 
nance has been maximized. In correspondence with this 
is the relative minimum of Tmoas around the 3-photon 
resonance, see Fig. IHlb). Interestingly enough, the time 
scale of the measurement time around this resonance is 
Tmoas ~ 10~^ X 27r/ri. Considering realistic numbers of 
a typical experimental set-up^, where fl is in the regime 
of a few GHz. we obtain a time scale of Tmoas ~ 100 ps 
for the nonlinear quantum detection scheme. This should 
be contrasted to the measurement time of Tmoas ~ 300 
ns obtained in Ref. 0. In between the multiphoton res- 
onances, the dependence of T^^cas on shows a rich 
structure including several singularities, which are sim- 
ply due to the several crossings of the two nonlinear re- 
sponse curves shown in Fig. [5]a), where Z?(a;ox) becomes 
zero, implying insufficient discrimination of the two qubit 
states. 

With this, we can evaluate the measurement efficiency, 
defined by the ratio Fmcas/r, with Fmeas = T;^Ls- This 
quantity sets the probability to infer the qubit state, 
based on the nonlinear response of the detector. We show 
the result for the efficiency of the measurement in Fig. 
[ni^c). Related to the multiphoton resonances in the de- 
tector, the efficiency also shows local maxima. For the 
discrimination power being optimized around the three- 
photon resonance, the measurement efficiency displays a 
clear local maximum [see Fig. IHc)]. Due to the small 
size of the relaxation rate F of the detector, the overall 
measurement efficiency is rather large in comparison to 
the detection set-up with a linear resonator,— ensuring 

rnicas/r 1. 



VII. EFFICIENCY OF THE MEASUREMENT 



VIII. CONCLUSIONS 



The measurement of the qubit state requires a coupling 
to the outer world, which clearly introduces noise to the 
qubit. In turn, the noisy detector yields measurement 
results, which are statistically distributed. This implies 
that several measurements have to be performed to ob- 
tain a reliable statistics. Hence, the relaxation time of the 
qubit state should not only exceed the typical relaxation 



To conclude, we have introduced a scheme for quantum 
state detection on the basis of a nonlinear detector which 
is operated in the regime of resonant few-photon transi- 
tions. Discrete multiphoton resonances in the detector 
can be used to infer the state of the parametrically cou- 
pled qubit via a state-dependent frequency shift of the 
detector's nonlinear response function. The multiphoton 
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FIG. 6. (Color online) (a) Relaxation rate F of the nonlinear 
quantum detector, (b) the measurement time Tmeas, and (c) 
the measurement efhciency Fmoas/r as a function of the ex- 
ternal frequency ojcx- The parameters are the same as in Fig. 
II 



resonances are well separated in the spectrum and sharp 
enough to allow for a good resolution of the qubit state. 

By analyzing key quantities of the detector, we have 
shown that the nonlinear few-photon detector can be op- 
erated efficiently, reliably, and with sufficiently weak back 
action. In fact, it can be efficiently tuned by tuning the 
amplitude of the ac bias current of the SQUID. Further- 
more, we have shown that the sharpness of the multi- 
photon resonances can be used to obtain an increased 
discrimination power as compared to the linear paramet- 
ric detection scheme. Clearly, the relaxation rate at a 
multiphoton resonance for the qubit becomes maximal, 
but in general remains very small. The measurement 
time around a multiphoton resonance can be tuned such 
that it becomes minimal. For realistic experimental pa- 
rameters, we find surprisingly small measurement times, 
allowing in principle for fast measurements. Moreover, 
the efficiency of the measurement, which takes the time 
to acquire enough information to infer the qubit state 
into account, also assumes large values, thus allowing for 
a reliable and highly efficient measurement of the qubit 
state. 

We have chosen realistic values for the involved model 
parameters such that an experimental realization of this 
quantum measurement scheme should become possible in 
the near future. The nonlinear detection scheme in the 
deep few-photon quantum regime offers thus the advan- 
tage of an increased discrimination power of more than 
98% (for our choice of realistic parameters), as com- 
pared to previous classical detection schemes based on 
the Joscphson bifurcation amplifier. 

A possible setup in order to realize the nonlinear few- 
photon detector could be the architecture used in a re- 
cent experiment*^ The low-temperature regime, where 
quantum noise effects are important, has already been 
reached. In order to operate in the regime of only few 
photons in the resonator, the sensitivity and stability of 
the devices might have still to be further increased. How- 
ever, no principle obstacles are apparent. 
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Appendix: Power spectrum 

To calculate the power spectrum of a driven out-off- 
equilibrium system is a non trivial task, since the time re- 
versal symmetry, which simplifies the calculation in equi- 
librium systems, is not given anymore. An elegant way 
to compute this in terms of correlation functions is pre- 
sented in Ref. [2^ The correlation function between two 
operators A and B is described as a mean value 

SAB{t,T)=TvseB{W{t + T)A{0)}. (A.l) 

Here, the trace is over the whole system-plus- 
bath, with the "density" operator 'W{t + t) = 
Uit + T,t)[BWit)]U^t + T,t) and J7(t + r,t) = 

exp{-iT jl^'' Htota\{t')dt'}, with iJtotai being the Hamil- 
tonian of the whole system-plus-bath, W the density op- 
erator of the total system, and T the time ordering op- 
erator. Furthermore, A and B arc in the Hcisenberg 
representation. 

In the regime of weak coupling to the environment, 
the reduced density operator Q{t -\- t) = TrB{yV(t + r)} 
evolves according to the master equation (|5]). 

In the superoperator notation, Va^.a'p' = — — 
£l3)Saa'Sl3l3' + J~-af3,a'i3' (LiouviUc supcropcrator) is rep- 
resented by a Af^ x Af^ supermatrix V, where Af is the 
number of effective states in the system. In the same way, 
the density operator formally is a Af^ dimensional 
column vector g{t). The solution of the master equation 
is reduced to an eigenvalue problem of the matrix V, as 

2?.v™ = F,„v™, vl-V^r,„vl, (A.2) 

where v™ and are the left and right eigenvectors, 
respectively, with eigenvalue Tm- 

In the superoperator notation, the master equation (|5]) 
is expressed as g(t) = V ■ g{t), and its solution is given 
by 

£»(t) = exp{2?t} • £»(t = io)- (A.3) 

In the regime of the RWA, and at low temperature the 
master equation ([8|) conserves the trace and the posi- 
tivity of the density operator, i.e., it assumes Lindblad 
form. Therefore, we can expand the solution of the mas- 
ter equation in terms of the right-eigenvector v™ 

Q{t)^Y.^"'c,nCxp{r„,t}, (A.4) 

m 

with Cm = vj„ • g{to)- Here, we have used the orthog- 
onality property v}^ ■ v™ = 6mm'- The corresponding 
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expression follows for the operator g{t + r), but with a 
different initial condition. It is easily understandable in 
the operator notation 

Qcpit + r) = <fi Bc'^'it) Q^,p,{t) e^"-^ 

(A.5) 

where Bapit) = {<t)a{t)\B\(t>0{t)) , and 

iym,,ap) is the 

operator representation of v™ (vm) in the quasienergy 
states. Considering an initial time in the stationary 
regime, i.e., gap{t) — > £>J^, and after averaging the initial 
time t over the period 27r/a;ox, the correlation function 
reads 



5'aB (t) ~ ^ Snm e 



Tjnt — inLJcxi 



(A.6) 



with 



al3 a'P'C 



where Aa/s^n and 



(A.7) 

are the coefBcients of the Fourier 



expansion of Aaf3{t) and Bap{t), respectively. The power 
spectrum is obtained directly from the Fourier transform 
ofEq. (1X61) . 
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